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ABSTRACT 

This article presents a new model for valuing financial contracts subject to credit risk and collateralization. 
Examples include the valuation of a credit default swap (CDS) contract that is affected by the trilateral 
credit risk of the buyer, seller and reference entity. We show that default dependency has a significant 
impact on asset pricing. In fact, correlated default risk is one of the most pervasive threats in financial 
markets. We also show that a fully collateralized CDS is not equivalent to a risk-free one. In other words, 
full collateralization cannot eliminate counterparty risk completely in the CDS market. 
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1 Introduction 


There are two primary types of models that attempt to describe default processes in the literature: 
structural models and reduced-form (or intensity) models. Many practitioners in the credit trading arena 
have tended to gravitate toward the reduced-from models given their mathematical tractability. 

Central to the reduced-form models is the assumption that multiple defaults are independent 
conditional on the state of the economy. In reality, however, the default of one party might affect the default 
probabilities of other parties. Collin-Dufresne et al. (2003) and Zhang and Jorion (2007) find that a major 
credit event at one firm is associated with significant increases in the credit spreads of other firms. Giesecke 
(2004), Das et al. (2006), and Lando and Nielsen (2010) find that a defaulting firm can weaken the firms in 
its network of business links. These findings have important implications for the management of credit risk 
portfolios, where default relationships need to be explicitly modeled. 

The main drawback of the conditionally independent assumption or the reduced-form models is 
that the range of default correlations that can be achieved is typically too low when compared with empirical 
default correlations (see Das et al. (2007)). The responses to correct this weakness can be generally 
classified into two categories: endogenous default relationship approaches and exogenous default 
relationship approaches. 

The endogenous approaches include the contagion (or infectious) models and frailty models. The 
frailty models (see Duffie et al. (2009), Koopman et al. (2011), etc) describe default clustering based on 
some unobservable explanatory variables. In valuations of contagion or infectious type models (see Davis 
and Lo (2001), Jarrow and Yu (2001), etc.), the assumption of conditional independence is relaxed and 
default intensities are made to depend on default events of other entities. Contagion and frailty models fill 
an important gap but at the cost of analytic tractability. They can be especially difficult to implement for 
large portfolios. 
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The exogenous approaches (see Li (2000), Laurent and Gregory (2005), Hull and White (2004), 
Brigo et al. (2011), etc) attempt to link marginal default probability distributions to the joint default 
probability distribution through some external functions. Due to their simplicity in use, practitioners lean 
toward the exogenous ones. 

Given a default model, one can value a risky derivative contract and compute credit value 
adjustment (CVA) that is a relatively new area of financial derivative modeling and trading. CVA is the 
expected loss arising from the default of a counterparty (see Brigo and Capponi (2008), Lipton and Sepp 
(2009), Pykhtin and Zhu (2006), Gregory (2009), Bielecki et al (2013) and Crepey (2015), etc.) 

Collateralization as one of the primary credit risk mitigation techniques becomes increasingly 
important and widespread in derivatives transactions. According the ISDA (2013), 73.7% of all OTC 
derivatives trades (cleared ad non-cleared) are subject to collateral agreements. For large firms, the figure 
is 80.7%. On an asset class basis, 83.0% of all CDS transactions and 79.2% of all fixed income transactions 
are collateralized. For large firms, the figures are 96.3% and 89.4%, respectively. Previous studies on 
collateralization include Johannes and Sundaresan (2007), Fuijii and Takahahsi (2012), Pitcrbarg (2010), 
Bielecki, et al (2013) and Hull and White (2014), Xiao (2017), etc. 

This paper presents a new framework for valuing defaultable financial contracts with or without 
collateral arrangements. The framework characterizes default dependencies exogenously, and models 
collateral processes directly based on the fundamental principals of collateral agreements. For brevity we 
focus on CDS contracts, but many of the points we make are equally applicable to other derivatives. CDS 
has trilateral credit risk, where three parties - buyer, seller and reference entity - are defaultable. 

In general, a CDS contract is used to transfer the credit risk of a reference entity from one party to 
another. The risk circularity that transfers one type of risk (reference credit risk) into another (counterparty 
credit risk) within the CDS market is a concern for financial stability. Some people claim that the CDS 
market has increased financial contagion or even propose an outright ban on these instruments. 

The standard CDS pricing model in the market assumes that there is no counterparty risk. Although 
this oversimplified model may be accepted in normal market conditions, its reliability in times of distress 


2 



has recently been questioned. In fact, counterparty risk has become one of the most dangerous threats to 
the CDS market. For some time now it has been realized that, in order to value a CDS properly, counterparty 
effects have to be taken into account (see ECB (2009)). 

We bring the concept of comvariance into the area of credit risk modeling to capture the statistical 
relationship among three or more random variables. Comvariance was first introduced to economics by 
Deardorff (1982), who used this measurement to correlate three factors in international trading. Furthermore, 
we define a new statistics, comrelation, as a scaled version of comvariance. Accounting for default 
correlations and comrelations becomes important in determining CDS premia, especially during the credit 
crisis. Our analysis shows that the effect of default dependencies on a CDS premium from large to small 
accordingly is i) the default correlation between the protection seller and the reference entity, ii) the default 
comrelation, iii) the default correlation between the protection buyer and the reference entity, and iv) the 
default correlation between the protection buyer and the protection seller. In particular, we find that the 
default comvariance/comrelation has substantial effects on the asset pricing and risk management, which 
have never been documented. 

There is a significant increase in the use of collateral for CDS after the recent financial crises. Many 
people believe that, if a CDS is fully collateralized, there is no risk of failure to pay. Collateral posting 
regimes are originally designed and utilized for bilateral risk products, e.g., interest rate swap (IRS), but 
there arc many reasons to be concerned about the success of collateral posting in offsetting the risk of CDS 
contracts. First, the value of CDS contracts tends to move very suddenly with big jumps, whereas the price 
movements of IRS contracts are far smoother and less volatile than CDS prices. Second, CDS spreads can 
widen very rapidly. Third, CDS contracts have many more risk factors than IRS contracts. In fact, our model 
shows that full collateralization cannot eliminate counterparty risk completely for a CDS contract. 

The rest of this paper is organized as follows: Pricing multilateral defaultable financial contract is 
elaborated on in Section 2; numerical results are provided in Section 3; the conclusions are presented in 
Section 4. All proofs and some detailed derivations are contained in the appendices. 
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2 Pricing Financial Contracts Subject to Multilateral Credit Risk 

We consider a filtered probability space (Q , J 7 , \ZP },, 0 , P ) satisfying the usual conditions, 
where Q denotes a sample space, J 7 denotes a a -algebra, P denotes a probability measure, and 
{jpf },> 0 denotes a filtration. 

In the reduced-form approach, the stopping (or default) time r, of firm i is modeled as a Cox arrival 
process (also known as a doubly stochastic Poisson process) whose first jump occurs at default and is 
defined by, 

r, = inf{/: \^h i (s,Z s )ds>H i \ (1) 

where h t (t) or h i (t,Z ! ) denotes the stochastic hazard rate or arrival intensity dependent on an exogenous 
common state Z,, and H i is a unit exponential random variable independent of Z f . 

It is well-known that the survival probability from time t to s in this framework is defined by 

Pi (t, s) := P t (T > s \ t > t,Z t ) = exp^- j h l (u)du) (2a) 

The default probability for the period (t, s ) in this framework is given by 

q t (t,s) := P t (T < s I t > t,Z t ) = l-pi(t,s) = l-expf- J h^ujdiA (2b) 

There is ample evidence that corporate defaults are correlated. The default of a firm’s counterparty 
might affect its own default probability. Thus, default correlation/dependence occurs due to the 
counterparty relations. 

The interest in the financial industry for the modeling and pricing of multilateral defaultable 
instruments arises mainly in two respects: in the management of credit risk at a portfolio level and in the 
valuation of credit derivatives. Central to the valuation and risk management of credit derivatives and risky 
portfolios is the problem of default relationship. 

Let us discuss a three-party case first. A CDS is a good example of a trilateral defaultable instrument 
where the three parties are counterparties A, B and reference entity C. In a standard CDS contract one party 
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purchases credit protection from another party, to cover the loss of the face value of a reference entity 
following a credit event. The protection buyer makes periodic payments to the seller until the maturity date 
or until a credit event occurs. A credit event usually requires a final accrual payment by the buyer and a 
loss protection payment by the protection seller. The protection payment is equal to the difference between 
par and the price of the cheapest to deliver (CTD) asset of the reference entity on the face value of the 
protection. 

A CDS is normally used to transfer the credit risk of a reference entity between two counterparties. 
The contract reduces the credit risk of the reference entity but gives rise to another form of risk: 
counterparty risk. Since the dealers arc highly concentrated within a small group, any of them may be too 
big to fail. The interconnected nature, with dealers being tied to each other through chains of OTC 
derivatives, results in increased contagion risk. Due to its concentration and interconnectedness, the CDS 
market seems to pose a systemic risk to financial market stability. In fact, the CDS is blamed for playing a 
pivotal role in the collapse of Lehman Brothers and the disintegration of AIG. 

For years, a widespread practice in the market has been to mark CDS to market without taking the 
counterparty risk into account. The realization that even the most prestigious investment banks could go 
bankrupt has shattered the foundation of the practice. It is wiser to face frankly the real complexities of 
pricing a CDS than to indulge in simplifications that have proved treacherous. For some time now it has 
been realized that, in order to value a CDS properly, counterparty effects have to be taken into account. 

Let A denote the protection buyer, B denote the protection seller and C denote the reference entity. 
The binomial default rule considers only two possible states: default or survival. Therefore, the default 
indicator y ■ for firm j (j = A or B or C) follows a Bernoulli distribution, which takes value 1 with default 

probability q ] , and value 0 with survival probability p . The marginal default distributions can be 

determined by the reduced-form models. The joint distributions of a multivariate Bernoulli variable can be 
easily obtained via the marginal distributions by introducing extra correlations. The joint probability 
representations of a trivariate Bernoulli distribution (see Teugels (1990)) are given by 
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Pooo -=P(Ya = ®’Y b = 0 ,Y C =0) = P a PbPc + Pc <j ab + Pb <j ac + Pa <j bc ~®abc (3 a ) 

Pioo := P(Y a =1 ,Y B =0 ,Y C =0) = q A P B Pc ~ Pc °ab ~ Pb^ac +c 1a <j bc + ®abc (3b) 

Pow := P(Y a = 0, Y b = 1, Y c = 0) = PaPbPc ~ Pc a ab + qB <J AC ~ Pa&bc + @abc (3c) 

Pool := P(Y a = QY b =0 ,Y C =1) = PaPbQc + Pc^ab ~PB a ac ~PA a Bc + @abc (3d) 

Pi 10 := P(Y A =b^c = ty = PaPbPc + Pc a AB ~PB a AC ~qA a BC ~ PABC (3 e ) 

Pioi -=P(Y a =1 ,Y b =0 ,Y C =1) = ^aPb^c ~Pc a ab + Pb <j ac ~qA a Bc — PABC (3f) 

Poll := P(Ya = 0,Y b =1 ,Y C = 1) = p aPbPc ~Pc a ab ~Pb (T ac + Pa^bc ~Pabc (3g) 

Pill := P(Ya = 1 'Y B =YY c =1) =qA<iB c lc + qc a AB + qB (7 AC + qA (7 BC + @ABC (3h) 

where 

e mc := E{(Y a - q A )(Y B -q B )(Y c -q c )) (3i) 


Equation (3) tells us that the joint probability distribution of three defaultable parties depends not 
only on the bivariate statistical relationships of all pair-wise combinations (e.g., cr,.) but also on the 

trivariate statistical relationship (e.g., d ABC ). d ABC was first defined by Deardorff (1982) as comvariance, 

who use it to correlate three random variables that are the value of commodity net imports/exports, factor 
intensity, and factor abundance in international trading. 

We introduce the concept of comvariance into credit risk modeling arena to exploit any statistical 
relationship among multiple random variables. Furthermore, we define a new statistic, comrelation, as a 
scaled version of comvariance (just like correlation is a scaled version of covariance) as follows: 
Definition 1: For three random variables X A , X B ,and X c ,let /i A , ju B ,and ju c denote the means of X A 
, X B ,and X c . The comrelation of X A , X B ,and X c is defined by 

r _ e[(X a - Ma )(Xb-Mb)(X c -M c )\ 

9 ABC - I ------ W 

P\X a -Ma\ xE \ x b -Mb\ xE\X c - ju c \ 

According to the Holder inequality, we have 
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I E((X a -fi A )(X R -H R )(X C - Mc ))\< E\(X a -ju a )(X b - Mb ){X c -M e)\ ( < 

<l]E\X A -m a \ 3 xE\X b -m b \' *E\X c -Met 

Obviously, the comrelation is in the range of [-1, 1], Given the com relation, Equation (3i) can be 


rewritten as 


0 ABC := E{(Y a - q A )(Y B - q B )(X c - q c )) = Cabc^E\Y a - q A | 3 x E\Y b - q B | 3 x E\Y C - 9c | 3 

r -----( 6 ) 

= £ abc^\P a 1 ! a(P a + c 1a)Pb c Ib(Pb + c 1b)Pc c 1c(Pc +c 1c)\ 

where E(Yf) = qj and e|f. -q^ = p q.(/>■ • q )), j=A, B, or C. 

If we have a series of n measurements of X A , X B , and X, written as x Ai , x Bi and x Ci where / = 
1,2,.. ,,n, the sample comrelation coefficient can be obtained as: 


CABC ~ 


_ Z'li (XAi - Ma )(■% - Mb )( x a - Be ) _ 

a/ZLKu x Z"=il x Bi -/^b | 3 x Z"=i \ x ct-Bc \ 3 


More generally, we define the comrelation in the context of n random valuables as 
Definition 2: For n random variables X lt X 21 ... t X n , let Mi denote the mean of X f where i=l,..,n. The 
comrelation of X lt X 1> ... > X n is defined as 


Cl2 ...« _ ' 


g[(X 1 -/q)(X 2 -/i 2 )--.(X,,-; 0 ] 
/£|X 1 -/i 1 |"x£|X 2 -/i 2 |"-x£|X„ -M„\" 


Correlation is just a specific case of comrelation where n = 2. Again, the comrelation is in 
the range of [-1, 1] according to the Holder inequality. 

2.1 Risky valuation without collateralization 

Recovery assumptions are important for pricing credit derivatives. If the reference entity under a 
CDS contract defaults, the best assumption, as pointed out by J. P. Morgan (1999), is that the recovered 
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value equals the recovery rate times the face value plus accrued interest 2 . In other words, the recovery of 
par value assumption is a better fit upon the default of the reference entity, whereas the recovery of market 
value assumption is a more suitable choice in the event of a counterparty default 3 . 

Let valuation date be t. Suppose that a CDS has m scheduled payments represented as 
Xj = -sNS(T i _ l ,T i ) with payment dates T x where S(T i l ,T i ) denotes the accrual factor for 

period , N denotes the notional/principal, and s denotes the CDS premium. Party A pays the 

premium/fee to party B if reference entity C does not default. In return, party B agrees to pay the protection 
amount to party A if reference entity C defaults before the maturity. We have the following proposition. 
Proposition 1: The value of the CDS is given by 




(9a) 


where t =T 0 and 


0(Tj, T j+1 ) l( y(T/4i )+x , +1 )> 0 ( K ( Tj, T j+l ) + l( y(r . +i)+x , +i )< 0 </> A ( Tj, T j+1 ) 


(9b) 


K (Tj , T ]+1 ) = {p A {Tj , T j+1 ) Pb {Tj , T J+1 )p c {Tj , T j+l ) + q A (T s , T J+1 )p B {Tj , T j+l )p c (Tj , T J+i )cp A (T j+1 ) 

+ Pa ( T j . T j + 1 )<7« (Tj,T J+l )p c (Tj,T j+ , )cp A (T J+1 ) + q A (Tj , T j+1 )q B (T,T J+l )p c (Tj , T j+l )(p AB (T j+1 ) 

+ Pc (Tj, Tj+i )a AB (Tj , T J+i )(l - cp A (T j+l ) - <p A (T j+1 ) + cp AB (T j+1 )) (9c) 

+ cr AC (Tj , T J+1 )[p B (Tj,T j+l )(l - (p A (T ) +l ))+ q B (T , T j+1 ){<p A (T j+1 ) - (p AB (T j+l ))] 

+ cr SC (Tj - T i+ 1 )[p A (Tj - T J+ 1 )(* “ Va ( T j+ i )) + d A ( T j . )(^a ( r y + i) “ C^+i ))] 

+ Oabc (Tj , r y+1 )(-1+^A (r ;+1 ) - ^ab (T j+i ) + ^A (T M ) )}d(Tj , r. +1 ) 


2 In the market, there is an average accrual premium assumption, i.e., the average accrued premium is half 
the full premium due to be paid at the end of the premium. 

3 Three different recovery models exist in the literature. The default payoff is either i) a fraction of par 
(Madan and Unal (1998)), ii) a fraction of an equivalent default-free bond (Jarrow and Turnbull (1995)), or 
iii) a fraction of market value (Duffie and Singleton (1999)). 
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K (Tj,T J+l ) = {p A (T j , T j+1 )p B (T, , T j+1 ) Pc (P,, T j+1 ) + q A (P,, T j+1 ) p B (T j ,P. +I )p c (T j , T j+l )<p B (T j+1 ) 

+ Pa (T'j > 1 )^b > ^ 7 + 1 ) Pc (Tj > Tj+i )<Pb (Tj+i ) + a ’ Tj+i ) c Ib (Tj , P,- +1 ) p c (P ; , 7V +1 )^ AB (Pj+i) 

+ Pc (T,, Tj+l)a AB (Tj ,P ,. +1 )(l - (p B (T j+ i ) - <p B (T j+l ) + cp AB (T j+I )) (9d) 

+ <^ c (Tj - r , +1 )[p B C Tj , P , +1 )(l ■- cp B (T j+i )) + q B C Tj , 7 . +1 )(p B (T m ) - (T J+1 ))] 

+ <?bc (Tj - T J+ i )[p A (Tj , T j+ i )(l - cp B (Tj )) + q A (Tj , Tj+i ){<p B (Tj+i ) - (p AB (Tj+i ))] 

+ 0(Tj , T j+l )(-1 + <p B (Tj+i ) - cp AB (Tj+i ) + cp B (Tj+i ))}D(Tj , T J+1 ) 

CXTj , Tj+i ) = {p A (Tj , T ]+ i )p B (Tj , Tj+i )q c (Tj , T j+l ) + q A (Tj , P ,. +1 )p B (7,, P, +I ) 9c (P,, P , +1 )p B (P, +I ) 

+ Pa (^j- > Tj+i )q B (Tj, Tj+i )q c (Tj, T j+i )cp B (T ; . +1 ) + q A (Tj, Tj +i )q B (Tj , T j+1 )q c (Tj, Tj +i )(p AB (T j +1 ) 

+ q c (Tj , Tj+i )a AB (T ), P J+I )(l - (p B (T J+1 ) - (p B (T J+1 ) + <p AB (T j+l )) (9e) 

— &ac (Tj>Tj+i)\p B (Tj , Tj +I )(l — (p B (Tj +l ))+ q B (Tj ,Tj +l )[(p B (T J+l ) — cp AB (P ;+ ,))] 

- (Tj ,P, + ,)[p A (P,, P , +1 )(l - cp B (Tj+l ))+q A (Tj , Tj+l )(p B (P, +1 ) - ^ AB (P, +1 ))] 

^ABC (Tj 5 )(l - <P R (Tj+i) + (p AB (Tj+1 ) - cp B (Tj+i))}D(Tj , T J+ i ) 

where P(P.,P. +1 ) = (iv(l-^ c (P. + 1 ))-a(P.,P. +1 )), «(P.,P, +1 ) = sNS(T s ,T)/2, and X, =-sNS(T r T j+l ). 


Proof: See the Appendix. 

We may think of 0(t,T ) as the risk-adjusted discount factor for the premium and Q(t,P) as the 
risk-adjusted discount factor for the default payment. Proposition 1 says that the pricing process of a 
multiple-payment instrument has a backward nature since there is no way of knowing which risk-adjusted 
discounting rate should be used without knowledge of the future value. Only on the maturity date, the value 
of an instrument and the decision strategy are clear. Therefore, the evaluation must be done in a backward 
fashion, working from the final payment date towards the present. This type of valuation process is referred 
to as backward induction. 

Proposition 1 provides a general form for pricing a CDS. Applying it to a particular situation in 
which we assume that counterparties A and B are default-free, i.e., Pj = 1, r/ ; = 0 , p kl = 0, and g ABC = 0 , 

where j-A or B and k, l-A, B, or C, we derive the following corollary. 

Corollary 1: If counterparties A and B are default-free, the value of the CDS is given by 


=z;:, £ ltn "o ^ 1-^ J + z z, ^itn ;: 2 o^ ^ J 

= El, E l D (T T ,) Pc (I T l )X i \f t ]+ x;, e[dU, T ) Pc (t, T , )q c (T ,, T)R(T, ,, T )\f t ] 


( 10 ) 


where ()(T ,,P) = D(T i t ,P) p c (T. ,,P.j; Q.(T i _ 1 ,T i ) = D(T i _ 1 ,T i )q c (T i _ 1 ,T i ). 
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The proof of this corollary becomes straightforward according to Proposition 1 by setting p kl =0, 

<p AB = o, W = 0 , Pj = 1, qj = o , Pc (t,T t ) = U^oP (T v T , + 0 , and D(t, 7j.) = • 

If we further assume that the discount factor and the default probability of the reference entity arc 
uncorrelated and the recovery rate cp c is constant, we have 

Corollary 2: Assume that i) counterparties A and B are default-free, ii) the discount factor and the default 
probability of the reference entity are uncorrelated; Hi) the recovery rate cp c is constant; the value of the 

CDS is given by 

V(t) = YZi P ( t ' T i)Pc(tJi-i)q c (T i _ 1 J i ){N(l-<p c )-a(T i _ l Jj)-Y!!l l P(tJ i )Pc(tJi)sNS(T i _ 1 JJ (11) 
where Pit.Tf = £ , [zXf,7’)|j^] denotes the bond price, p c it, 7}) = E\p c (t,T t ) \j r t ] , q c {t,T i ) = \-p c (t,T i ), 
Pit, T t _ } )q(T^ ,f) = p(t , T t _ x ) - pit, T i ) . 

This corollary is easily proved according to Corollary 1 by setting t[a') / | j^] = /J[X | ]/f [k| ] when 

X and Y arc uncorrelated. Corollary 2 is the formula for pricing CDS in the market. 

Our methodology can be extended to the cases where the number of parties n > 4. A generating 
function for the (probability) joint distribution (see details in Teugels (1990)) of 72 -variate Bernoulli can be 
expressed as 



where 0 denotes the Kronecker product; p in) = \p ( f ] } and cj in) = \ are vectors containing 2 


components: pf = p ki ^... K , * = l+2,"=i*/ 2< ' ’ k i G i 0 ’ 1 }’ =< 




2.2 Risky valuation with collateralization 

Collateralization is the most important and widely used technique in practice to mitigate credit risk. 
The posting of collateral is regulated by the Credit Support Annex (CSA) that specifies a variety of terms 
including the threshold, the independent amount, and the minimum transfer amount (MTA), etc. The 
threshold is the unsecured credit exposure that a party is willing to bear. The minimum transfer amount is 
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the smallest amount of collateral that can be transferred. The independent amount plays the same role as 
the initial margin (or haircuts). 

In a typical collateral procedure, a financial instrument is periodically marked-to-market and the 
collateral is adjusted to reflect changes in value. The collateral is called as soon as the mark-to-market 
(MTM) value rises above the given collateral threshold, or more precisely, above the threshold amount plus 
the minimum transfer amount. Thus, the collateral amount posted at time t is given by 


C(t) 


(V(t)-H(t) 

to 


if V(t) > Hit) 
otherwise 


(13) 


where H{t) is the collateral threshold. In particular, H(t) = 0 corresponds to full-collateralization 4 ; H >0 


represents partial/under-collateralization; and H< 0 is associated with over-collateralization. Full 


collateralization becomes increasingly popular at the transaction level. In this paper, we focus on full 
collateralization only, i.e., C(t) =V(t). 


The main role of collateral should be viewed as an improved recovery in the event of a counterparty 
default. According to Bankruptcy law, if there has been no default, the collateral is returned to the collateral 
giver by the collateral taker. If a default occurs, the collateral taker possesses the collateral. In other words, 
collateral does not affect the survival payment; instead, it takes effect on the default payment only. 

According to the ISDA (2013), almost all CDSs are fully collateralized. Many people believe that 
full collateralization can eliminate counterparty risk completely for CDS. 

Collateral posting regimes are originally designed and utilized for bilateral risk products, e.g., IRS, 
but there are many reasons to be concerned about the success of collateral posting in offsetting the risks of 
CDS contracts. First, the values of CDS contracts tend to move very suddenly with big jumps, whereas the 
price movements of IRS contracts are far smoother and less volatile than CDS prices. Second, CDS spreads 


4 There are three types of collateralization: Full-collateralization is a process where the posting of collateral 
is equal to the current MTM value. Partial/under-collateralization is a process where the posting of collateral 
is less than the current MTM value. Over-collateralization is a process where the posting of collateral is 
greater than the current MTM value. 
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can widen very rapidly. The amount of collateral that one party is required to provide at short notice may, 
in some cases, be close to the notional amount of the CDS and may therefore exceed that party’s short-term 
liquidity capacity, thereby triggering a liquidity crisis. Third, CDS contracts have many more risk factors 
than IRS contracts. 

We assume that a CDS is fully collateralized, i.e., the posting of collateral is equal to the amount 
of the current MTM value: C(t) = V(t). For a discrete one-period (t, u ) economy, there are several possible 

states at time u: i) A, B, and C survive with probability p 000 . The instrument value is equal to the market 
value V(u) ; ii) A and B survive, but C defaults with probability p 00l . The instrument value is the default 
payment R(u ); iii) For the remaining cases, either or both counterparties A and B default. The instrument 
value is the future value of the collateral V(t)/ D(t,u) (Flere we consider the time value of money only). 
The value of the collateralized instrument at time t is the discounted expectation of all the payoffs and is 
given by 

V ( t ) = E{D(t, u)[p 000 (t, ujV (u) + p 001 (t, u)R(u) 

+ ( p l00 (t, u) + p 0l0 (t, u ) + p x 10 (t, u) + p l0l (t, u) + p 0l , {t, u) + p ni it, u )) V (, t) / D(t, u)]\f t } (14a) 

- E{ [D(t , u)(p m) (t, u)V (s) + Pool (, t, u)R(u))+ (l - Poooit, u ) - p 0(n (t, uj)V (f)] \f t } 

or 

^[(^ > a(1’ m )/ , b(Tm) + a AB (t, *0)1 fMt) 

= E§D(t, u)(p c (t, u)V (m) + q c (t, u)R(u))(p A (t, u)p B (t, u ) + a AB (t, u )) (14b) 

+ D(t,u)(v (m) — R(u)){p B (t , u)a AC (t , u) + p A (t ,u)cr BC ( t,u ) — 0 ABC (t,u 

If we assume that (p A (t,u)p B (t,u) + and D(t,u){p c (t,u)V(u) + q c (t,u)R(u)) are 

uncorrelated, we have 

V(t) = V F (t) +% ABC (t,u)/tf/ AB (t,u) (15a) 

where 

V F (t) = E{D(t,u)[p c (t,uW(u) + q c (t,u)R(u)]\jF t } (15b) 

^ AB (t ’ U ) = E\[p A (t,u)p B (t,li) + <7 AB it,u)]\jF t } (15c) 
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^abc^Tu) ~ E{\D(t,u)(p B (t,u)(TAc(t’ ,A ) + Pa (0 u ) cr nc (0 M ) @abc (U))(l/(«)-%)|/J (15d) 

The first term (f) in equation (15) is the counterparty-risk-free value of the CDS and the second 
term is the exposure left over under full collateralization, which can be substantial. 

Proposition 2: If a CDS is fully collateralized, the risky value of the CDS is NOT equal to the counterparty- 
risk-free value, as shown in equation (15). 

Proposition 2 or equation (15) provides a theoretical explanation for the failure of full 
collateralization in the CDS market. It tells us that under full collateralization the risky value is in general 
not equal to the counterparty-risk-free value except in one of the following situations: i) the market value 
is equal to the default payment, i.e., V(u) = R(u) ; ii) firms A, B, and C have independent credit risks, i.e., 

Pi! =0 and g ABC = 0 ; or iii) the following relationship holds p B cr AC + p A cr BC = 0 ABC . 

3 Numerical Results 

Our goal in this section is to study the quantitative relationship between CDS premia and the credit 
quality of counterparties and reference entities, including the default correlations and comrelations. 

In our study, we choose a new 5-year CDS with a quarterly payment frequency. Two counterparties 
are denoted as A and B. Counterparty A buys a protection from counterparty B. All calculations arc from 
the perspective of party A. By definition, a breakeven CDS spread is a premium that makes the market value 
of a given CDS at inception zero. 

The current (spot) market data are shown in Table 1 provided by FinPricing (2013). Assume that 
the reference entity C has an “A+200bps” credit quality throughout this subsection. The 5-year 
counterparty-risk-free CDS premium is 0.027 (equals the 5-year ‘A’ rated CDS spread in Table 1 plus 200 
basis points). 

Since the payoffs of a CDS arc mainly determined by credit events, we need to characterize the 
evolution of the hazard rates. Flere we choose the Cox-Ingersoll-Ross (CIR) model. The CIR process has 
been widely used in the literature of credit risk and is given by 
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dh t = a(b - h t )dt + crJh^dW l 


( 16 ) 


where a denotes the mean reversion speed, b denotes the long-term mean, and cr denotes the volatility. 


Table 1: Current/spot market data 

This table displays the current (spot) market data used for all calculations in this paper, including the term 
structure of continuously compounded interest rates, the term structure of A-rated breakeven CDS spreads, 


and the curve of at-the-money caplet volatilities. 


Term (days) 

31 

91 

182 

365 

548 

730 

1095 

1825 

2555 

3650 

5475 


0.002 

0.002 

0.002 

0.004 

0.007 

0.010 


0.024 

0.030 

0.035 

0.040 

Interest Rate 

8 

7 

9 

3 

1 

2 

0.016 

9 

6 

5 

5 

Credit Spread 

0.004 

0.004 

0.004 

0.004 

0.004 

0.005 

0.005 

0.007 

0.007 

0.009 

0.010 

2 

2 

2 

5 

9 

2 

8 


9 

1 

6 

Caplet 

0.326 

0.331 

0.337 

0.350 

0.364 

0.377 

0.308 

0.247 

0.214 

0.167 

0.163 

Volatility 

7 


6 

9 

1 

3 


3 

1 

8 

4 


Table 2: Risk-neutral parameters for CIR model 


This table presents the risk-neutral parameters that are calibrated to the current market shown in Table 1. 
‘A+lOObps’ represents a ‘100 basis points’ parallel shift in the A-rated CDS spreads. 


Credit Quality 

A 

A+100bps 

A+200bps 

A+300bps 

Long-Term Mean cl 

0.035 

0.056 

0.077 

0.099 

Mean Reverting Speed b 

0.14 

0.18 

0.25 

0.36 

Volatility cr 

0.022 

0.028 

0.039 

0.056 


The calibrated parameters arc shown in table 2. We assume that interest rates arc deterministic and 
select the regression-based Monte-Carlo simulation (see Longstaff and Schwartz (2001)) to perform risky 
valuation. 
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We first assume that counterparties A, B , and reference entity C have independent default risks, i.e., 
Pab = Pac = Pbc = ( Pab - Cabc = 0, an d examine the following cases: i) B is risk-free and A is risky; and ii) 
A is risk-free and B is risky. We simulate the hazard rates using the CIR model and then determine the 
appropriate discount factors according to Proposition 1. Finally we calculate the prices via the regression- 
based Monte-Carlo method. The results are shown in Table 3 and 4. 


Table 3: Impact of the credit quality of the protection buyer on CDS premia 

This table shows how the CDS premium increases as the credit quality of party A decreases. The 1st data 
column represents the counterparty-risk-free results. For the remaining columns, we assume that party B is 
risk-free and party A is risky. ‘A+100bps’ represents a ‘100 basis points’ parallel shift in the A-rated CDS 
spreads. The results in the row ‘Difference from Risk-Free’ = current CDS premium - counterparty-risk- 


free CDS premium. 


Credit Quality 

Party A 

~ 

A 

A+100bps 

A+200bps 

A+300bps 

Party B 

~ 

" 

" 

" 


CDS premium 

0.027 

0.02703 

0.02708 

0.02713 

0.02717 

Difference from Risk-Free 

0 

0.003% 

0.008% 

0.013% 

0.017% 


Table 4: Impact of the credit quality of the protection seller on CDS premia 

This table shows the decrease in the CDS premium with the credit quality of party B. The 1st data column 
represents the counterparty-risk-free results. For the remaining columns, we assume that party A is risk-free 
and party B is risky. ‘A+100bps’ represents a ‘100 basis points’ parallel shift in the A-rated CDS spreads. 
The results in the row ‘Difference from Risk-Free’ = current CDS premium - counterparty-risk-free CDS 


premium. 


Credit Quality 

Party A 

Party B 

A A+100bps A+200bps A+300bps 

CDS premium 

Difference from Risk-Free 

0.027 0.02695 0.02687 0.0268 0.02672 

0.00% -0.005% -0.013% -0.020% -0.028% 
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From table 3 and 4, we find that a credit spread of about 100 basis points maps into a CDS premium 
of about 0.4 basis points for counterparty A and about -0.7 basis points for counterparty B. The credit impact 
on the CDS premia is approximately linear. As would be expected, i) the dealer’s credit quality has a larger 
impact on CDS premia than the investor’s credit quality; ii) the higher the investor’s credit risk, the higher 
the premium that the dealer charges; iii) the higher the dealer’s credit risk, the lower the premium that the 
dealer asks. Without considering default correlations and comrelations, we find that, in general, the impact 
of counterparty risk on CDS premia is relatively small. This is in line with the empirical findings of Arora, 
Gandhi, and Longstaff (2009). 


Figure 1: Impact of default correlations and comrelation on CDS premia 

Each curve in this figure illustrates how CDS premium changes as default correlations and comrelation 
move from -1 to 1. For instance, the curve ‘cor BC’ represents the sensitivity of the CDS premium to 
changes in the correlation p BC when p AB = p AC = C ABC = 0. 


Impact of Default Correlations and Comrelation on CDS Premia 



correlation or comrelation 


-♦— cor_AB 
cor_AC 
-k — cor_BC 
■■— comr ABC 
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Next, we study the sensitivity of CDS premia to changes in the joint credit quality of associated 
parties. Sensitivity analysis is a very popular way in finance to find out how the value and risk of an 
instrument/portfolio changes if risk factors change. One of the simplest and most common approaches 
involves changing one factor at a time to see what effect this produces on the output. We are going to 
examine the impacts of the default correlations p AB , p AC , p BC , and the comrelation C' Afsc separately. 
Assume that party A has an ‘A+100bps’ credit quality and party B has an ‘A’ credit quality. The 5-year 
risky CDS premium is calculated as 0.02703. 

Assume <p AB =0.5. The impact diagrams of the default correlations and comrelation are shown in 
Figure 1. From this graph, we can draw the following conclusions: First, the CDS premium and the default 
co rre I a t i o n s/c o m relation have a negative relation. Intuitively, a protection seller who is positively correlated 
with the reference entity (a wrong way risk) should charge a lower premium for selling credit protection. 
Next, the impacts of the default correlations and comrelation are approximately linear. Finally, the 
sensitivity slopes of the CDS premium to the default correlations and comrelation are -0.06 to p AB ; -0.09 
to p AC ; -53 to p BC ; and -14 to £ AHC . Slope measures the rate of change in the premium as a result of a 
change in the default dependence. For instance, a slope of -53 implies that the CDS premium would have 
to decrease by 53 basis points when a default correlation/comrelation changes from 0 to 1. 

As the absolute value of the slope increases, so does the sensitivity. The results illustrate that p BC 
has the largest effect on CDS premia. The second biggest one is £ AHC ■ The impacts of p AB and p AC are 
very small. In particular, the effect of the comrelation is substantial and has never been studies before. A 
natural intuition to have on CDS is that the party buying default protection should worry about the default 
correlations and comrelation. 

4 Conclusion 

This article presents a new valuation framework for pricing financial instruments subject to credit 
risk. In particular, we focus on modeling default relationships. 
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To capture the default relationships among more than two defaultable entities, we introduce a new 
statistic: comrelation, an analogue to correlation for multiple variables, to exploit any multivariate statistical 
relationship. Our research shows that accounting for default correlations and comrelations becomes 
important, especially under market stress. The existing valuation models in the credit derivatives market, 
which take into account only pair-wise default correlations, may underestimate credit risk and may be 
inappropriate. 

We study the sensitivity of the price of a defaultable instrument to changes in the joint credit quality 
of the parties. For instance, our analysis shows that the effect of default dependence on CDS premia from 
large to small is the correlation between the protection seller and the reference entity, the comrelation, the 
correlation between the protection buyer and the reference entity, and the correlation between the protection 
buyer and the protection seller. 

The model shows that a fully collateralized CDS is not equivalent to a risk-free one. Therefore, we 
conclude that collateralization designed to mitigate counterparty risk works well for financial instruments 
subject to bilateral credit risk, but fails for ones subject to multilateral credit risk. 

Appendix 

Proof of Proposition 1. Let t = T 0 . On the first payment date 7j, let V (7j) denote the market 
value of the CDS excluding the current cash flow X, . There are a total of eight (2 3 = 8) possible states 
shown in Table A_l. 


Table A_l. Payoffs of a trilaterally defaultable CDS 

This table shows all possible payoffs at time 7j. In the case of V(T 1 ) + X l > 0 where V (7j) is the market 

value excluding the current cash flow X i , there are a total of eight (2 3 = 8) possible states: i) A, B, and C 
survive with probability p 000 . The instrument value equals the market value: P(7j) + X 1 . ii) A defaults, but 
B and C survive with probability p 100 . The instrument value is a fraction of the market value: 
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(p li (T l )(y(T t ) + A',) where <p B represents the non-default recovery rate of party B 5 . <p B = 0 represents the 
one-way settlement rule, while <p B =1 represents the two-way settlement rule, iii) A and C survive, but B 
defaults with probability p 0 io • The instrument value is given by (p,AT l )(V(7j) + X,) where (p B represents 
the default recovery rate of defaulting party B. iv) A and B survive, but C defaults with probability p 0 oi • 
The instrument value is the default payment: R(T 0 ,T t ) ■ v) A and B default, but C survives with probability 
Pi I,,. The insti ument value is given by (p AK (T i )(VT7j) + A,) where <p AB denotes the joint recovery rate 
when both parties A and B default simultaneously, vi) A and C default, but B survives with probability p lQl 
. The instrument value is a fraction of the default payment: <p R (T)R(T (j ,T l ). vii) B and C default, but A 
survives with probability p on , The instrument value is given by ip,,(T)A > (T il ,T 1 ). viii) A, B , and C default 
with probability p m . The instrument value is given by <p AP (T)R(T [t ,J \). A similar logic applies to the case 
of V(r 1 ) + X 1 <0. 


Status 

Probability 

Payoff if V(T 1 ) + X l >0 

Payoff if V(T 1 ) + X 1 <0 

o 

II 

o 

II 

o 

II 

>■? 

P 000 

V(T t ) + X t 

V(T t ) + X t 

o 

II 

O" 

II 

1 ^ 

II 

Pi 00 

P b (T i )(V(T i ) + a,) 

( P.\ (T\ )(V(T \) + A,) 

o 

II 

II 

o' 

II 

Poio 

VsKivcrj + xJ 

«p,(7j)(V^) + A,) 

II 

o' 

II 

o' 

II 

Pool 

R(T () ,T t ) 

R(T 0 , 7j) 

Y a =1,Y b =1,Y c =0 

Pi 10 

vMivcro + x,) 

Va B (T,){V{T x ) + X x ) 

ii 

to 

II 

o 

II 

Pm 

<p B {T)R{T 0 , 7j) 

<P r (T)R{T () ,T ] ) 

II 

II 

o' 

II 

Pon 

<Pb(T)R{T, o.t;) 

cp B (T)R(T n ,T l ) 


5 There are two default settlement rules in the market. The one-way payment rule was specified by the 
early ISDA master agreement. The non-defaulting party is not obligated to compensate the defaulting party 
if the remaining market value of the instrument is positive for the defaulting party. The two-way payment 
rule is based on current ISDA documentation. The non-defaulting party will pay the full market value of 
the instrument to the defaulting party if the contract has positive value to the defaulting party. 
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Y a =1,Y b = l,Y c =1 


Pi i 


9 AK (T)R{T {) ,T,) 


( Pab(T)R(T u ,T 1 ) 


The risky price is the discounted expectation of the payoffs and is given by 

v(0 = e{[i (y(r 1 )x 1 >o((f’ooo (T’o ,T x ) + p m (T 0 ,T l )p B (T 1 ) + p 010 (P 0 ,T 1 )<p„(T 1 ) +p no (T 0 , T x )(p AB (T x ))(V (^) + X, )| ) 

+ V(7i)+x,)<o ((PoooCTo^i) + p m (T 0 ,T 1 )<p A (T l ) + p mo (T 0 ,T 1 )<p A (T l ) + p no (T 0 ,T 1 )<p AB (T l )Xv(T 1 ) + X x )\f t ) 

+ (P00l( T 0’ T l)+ Pl0l( T 0’ T lWs( T l) + Pou( T 0’ T i)9b( T i)+ Pni( T 0’ T l)9AB( T l)) R ( T s’ T l)\ft] D (^ ,r )} 

= e{[O(T 0 ,T l )(V(T 1 ) + X l ) + Q(P 0 ,T ] )R(T , T x )]\f t } 


(Ala) 


where 


O(T 0 ,T]} — l(v )+Xj)>o Qb (To > 7]) + l ( y ( Ti )+Xl ><o Qa (To > 7)) 


(Alb) 


K(To’ T i) = {Pa(To >T'i )Pa(T’o>7) )Pc (T’o> T i ) + c l A ( T o »7))Ps(T'o>(T’o*7 )Va(7)) 

+ Pa (r 0 ,7) )? B (T’o, 7) )p c (T 0 ,7; )p A (7)) + q A (T 0 ,7) )c/ B (7^, 7) )p c (7^, 7) )^ AS (7)) 

+ Pc (To .7) )cr AB (7),, 7) )(l - y> A (7)) - (P,) + <p AB (7))) (Ale) 

+ Cac (T’o, T , )[p fi (P 0 ,P, )(l - ^ A (P,)) + q B (T 0 ,7) )(^ A (P,) - (p AB (7' 1 ))] 

+ Osc(T’o,7) )[p A (T 0 ,7))(l - cp A (7j)) + q A (7' 0 ,7' 1 )(p A (7)) - (p AB (P,))] 

+ R ABC (T 0 ,7j)(-1 + cp A (7j) - ^ AB (7)) + ^ A (7) ))\D(T 0 ,7' 1 ) 

0b( T o’ T i ) = {Pa(T’ o>7) )p s (P 0 ,P,)p c (P 0 , T y ) + q A (T 0 ,T x )p B (T 0 ,T l )p c (T 0 ,T x )(p B (T x ) 

+ p A (T x , T x )q B (T q , T x )p c (T 0 ,P, )(p B (P,) + q A (T 0 ,P, Vy B (P 0 ,T { )p c (T 0 ,T x )cp AB (T t ) 

+ Pc (T’o >7) )cr AB (P 0 ,P,)(l-p B (7)) - <p B (T x ) + ^ AB (Pj)) (Aid) 

+ ct ac (P 0 ,Pj )[/> a (P 0 , P, )(l ~V B (P,)) ■+ q B (T 0 ,P, )(<p B (P,) - ^ AB (P,))] 

+ ct sc (P 0 , T { )[p A (P 0 ,P, )(l - cp B (T x ))+q A (T 0 ,P, )fe(P,) - ^ AS (P,))] 

+ #asc (7’o»7’ 1 )(— 1 + ^ b (Pi)-^ab(7’i) + ^b(7’ 1 ))}P>(P 0 ,P 1 ) 

Q(P 0 , P) = {p A (P 0 , T)p B (T 0 , T)q c (P 0 , P) + ^ A (P 0 , T)p B (P 0 , P)^ c (P 0 , P)^ B (P,) 

+ Pa (TofT)q B (T 0 ,T)q c (T 0 ,T)<p B (T x ) + q A (T 0 ,T)q B (T 0 ,T)q c (T 0 ,T)cp AB (T x ) 

+ q c (T 0 , T)ct ab (P 0 , P)(l - <p B (P,) - (P,) + <3 AB (P,)) (Ale) 

- <r ac (To . 7’)[p a (T’o, 7')(l -^«(7’ l ))+ (P 0 , P)(^ B (Pi) - (p AB (T x ))] 

- a BC (P 0 , T)[p A (T q , P)(l -<p B (T x ))+ q A (P 0 , T){(p B (P,) - (p AB (P,))] 

+ ^ ASC (P 0 , P)(l - B (P,) + (T’o) - (T’o ))}P>(P 0 , P) 

Similarly, we have 

y(p 1 ) = p{[o(p 1 ,p 2 )(x 2 +y(p 2 ))+Q(p 1 ,p 2 )P(p 1 ,p 2 )]|j r J (A2) 
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Note that O(T 0 ,T t ) is J r T] -measurable. By definition, an J r T] -measurable random variable is a 


random variable whose value is known at time 7). According to taking out what is known and tower 

properties of conditional expectation, we have 

V{t) = E{[O(T 0 ,T l )( X i + V(T l ))+kl(T 0 ,T l )R(T Q ,T l )]\f t }=E[O(T 0 ,T l )X l \f t ] 

+ E[kl(T 0 , T { )R(T 0 , )|j f ] + EpiTgJ, )E^[0(T l ,T 2 ){X 2 +V(T 2 ))+ Q(F,, T 2 )R(T ] , T 2 )] } | J t } 

-ZL, £ [(nr>(r.r.,)) x J+lL £ [(n"„ 0 <r.r..)) :i <r-..r)«<r-..r). 

+£ lnU 0 (r.r..)V(r)|j,. 


By recursively deriving from T 2 forward over T m , where V(T m ) = X m , we have 

v w=S‘, £ [(n;><r.r.,)) x ^J+Sr., £ [tnr> (7 'rr.r)h(r-,.r)*(r-,.r)|>; 


(A4) 
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